IMPERIAL/TP/95-96/56 



Nonperturbative calculation of symmetry breaking in quantum 

field theory 

Carl M. Benderffl 

Blackett Laboratory, Imperial College, London SW7 2BZ, United Kingdom 

Kimball A. Miltoni 

Department of Physics and Astronomy, University of Oklahoma, Norman, OK 73019, USA 

(February 1, 2008) 

Abstract 

A new version of the delta expansion is presented, which, unhke the con- 
ventional delta expansion, can be used to do nonperturbative calculations in 
a self-interacting scalar quantum field theory having broken symmetry. We 
calculate the expectation value of the scalar field to first order in 5, where 5 
is a measure of the degree of nonlinearity in the interaction term. 
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In a private discussion with D. Bessis we learned about the remarkable non-Hermitian 
Hamiltonian 

H = p^/2 + ix^. (1) 

The eigenvalues of this Hamiltonian are all real and positive! (We have no proof of this 
property, but there is ample numerical evidence.) This Hamiltonian is intriguing because it 
suggests a novel way to apply the delta expansion to quantum field theories having a broken 
symmetry. In this paper we introduce our new version of the delta expansion and use it to 
calculate the (nonzero) value of (0) in a scalar quantum field theory [see Eq. (PT|)]. 

The delta expansion is a Taylor series in powers of 5, where 5 measures the degree of 
nonlinearity of an interaction term For a scalar quantum field theory, the conventional 
approach has been to introduce the parameter 5 into the Euclidean Lagrangian density 
by 

C = (V0)V2 + m20V2 + g{<p'f+'. (2) 

The advantage of the delta expansion is that it is nonperturhative in the coupling constant 
g and mass m and has a nonzero radius of convergence. The conventional delta expansion 
has been used to study renormalization [|], supersymmetry [H, local gauge invariance p, 
stochastic quantization 0, and finite-temperature field theory 0. 

A drawback of the Lagrangian density in Eq. (H) is that it becomes a ^fj^l^ theory rather 
than a gclr' theory when 5 = \. The g\(l)\^ theory is symmetric under —(p and cannot 
exhibit symmetry breaking. In general, the conventional delta expansion is unsuitable for 
studying theories with symmetry breaking. 

In the past it was assumed that to satisfy the physical requirement that the Hamiltonian 
be bounded below for all 6, the parameter 6 should appear in the Lagrangian density (0) 
as the exponent of 0^ and not of 0. This assumption appears to be even more reasonable if 
one expands the interaction term in Eq. (Q) in powers of 6: 

9i<p'Y^'=9^'f:^-^Mr- (3) 

A;=0 
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One would think that the argument of the logarithm in Eq. (j^) should be positive to avoid 
complex numbers. 

However, H in Eq. ([1|) describes a theory having a positive spectrum. This suggests a 
surprising new version of the delta expansion in which we replace Eq. by 

£ = (V0) 72 + mVV2 - g{i(t)f+^. (4) 

As with Eq. (j^), 5 = gives free field theory. At 5 = 1 we recover (when the dimension of 
space-time is one and the bare mass vanishes) the theory described by Eq. (|l]). 

We believe that for Re5 > —2 the quantum-mechanical theory has a positive real 
spectrum. Here, we verify this to first order in 5. The shift in the nth energy level of 
H = p^/2 — {ixY~^^ is the expectation value of 6x'^\og{ix) in the unperturbed harmonic- 
oscillator basis \n): 

AEn = {n\5x'^ log{ix)\n) 

/oo 
dx ^Ij'^{x)x'^ \og{ix) 
-oo 

= 5 + J ^ dx iIj'^{x)x'^ \og{ix) 



oo 



S / dxip^{x)x [log(— ix) + log(ix)] 

26 / log(x) 
Jo 

S2~^/^[ar, - {2n + 1)7 - {5n + 5/2) log 2]. (5) 



[n/2+1/2] ^ 



The sequence a„ = 2, 8, 16, 74/3, ... is given by 

a„ = 4n + 4 + 4[n/2] + (4n + 2) E o/. 1 ' 

k=o - J- 

where [a] is the greatest integer in a. (In the above calculation we choose the branch of 
the logarithm to lie along the negative axis and choose the sheet of the Riemann surface for 
which logl = 0.) Although the Hamiltonian is not self-adjoint, the 0(5) shift in the energy 
levels is reall This property persists to all orders in 5 [Q. 

We will now show that in any Euclidean dimension and for all 6 > —2 (apart from a 
special discrete set of negative values of 6) , Eq. (^) exhibits symmetry breaking. This result 



is especially interesting for the case 5 = 2 because, although the resulting cj)^ theory appears 
symmetric under (f) —>■ —(f), symmetry breaking persists] 

We begin by considering a massless zero- dimensional scalar quantum field theory, whose 
vacuum-persistence amplitude Z and A^-point Green's function Gn are 

Z = r rfxe(-)^"^ Gj, = U^ dxe^-^'^'x^. (6) 

J —CO Zj J —oo 

To evaluate Gat we break the integral into two pieces: 



/U /"OO 
dxe^'^^'^'x"" + / dx e^'^^'^' x"" 
-oo Jo 



2(^:^^fj''Tl] rdxe^-^'^'x-. (7) 







^ilm (if N odd) 

This integral exists if —1 < Re 5 < 1. We let x^~^^ = s and use standard Fresnel contour 
integral methods to get 
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zGn = -i-^r / r . (8) 

^ ^ ' N +1 / \ 2 + 5 J ^ ^ 

Using this result we calculate Gn- 

The nonvanishing of G2N+1 in E<1- (|^) confirms that there is symmetry breaking. 

Let us examine the case 5 = 2 (0^ theory) in depth. Since the one-point Green's func- 
tion Gi for this model is nonvanishing, the symmetry remains permanently broken as 5 
approaches 2. The Green's functions for this theory are Gn = (— 'i)^r(3/4)/r(3/4 — N/A). 
In the presence of an external source J, the vacuum functional Z[J] for the theory has a 
convergent Taylor series in J: 

00 jNr^ -TT °° ( 



N=0 



2^oiV!r(^)- 



Note that Z[J] obeys the third-order differential equation 

Z'"[J] + JZ[J]/A = 0. (11) 



This is the D = form of the functional Schwinger- Dyson equation for the theory. If 
J is rotated by e*'^/^, Eq. ( p^T]) becomes the differential equation satisfied by the vacuum 
functional Z[J] for a conventional zero-dimensional 0^ theory in the presence of an external 
source: 



Z[J] = dxe-^'+-^\ (12) 

oo 



This model does not exhibit symmetry breaking. 

The Schwinger-Dyson equation ([TTD does not specify the vacuum functional uniquely 
PJTOH. Indeed, there are four distinct solutions to the rotated version of Eq. ([TT|): 



Z,[J]= / dxe-^^'^ (j = l, 2, 3, 4). (13) 

JCj 

Here, Cj are contours in the complex-x plane; Ci joins ioo to +oo, C2 joins +00 to — ioo, 
C3 joins —ioo to —00, and C4 joins ~oo to zoo. (The four solutions Zj[J] are not linearly 
independent; their sum is zero.) The vacuum functional Z[J] in Eq. (|T2]) is given by Z^[J] + 
Zi[J]. The vacuum functional for the symmetry-broken theory in Eq. (|IU|) is e*'^/^Z'i[e'^*'^/^ J]. 
The symmetry-broken theory corresponding to 5 = 2 in Eq. (|^) is distinguished from the 
conventional unbroken theory in Eq. ([T^) in that the sign of the term is reversed (reflecting 
the factor of e'^*'^/^ in the source term above). 

Next, let us examine the analytic behavior of Eq. (P) as a function of complex 5. The delta 
expansion of the one-point Green's function Gi, whose nonvanishing signals the presence of 
symmetry breaking, is 

Gi = -«v^((5/2 + (5^(7-2-21og2)/8-5^[48 + 7r2 

-6(7-4-2 log 2)2]/384 + ...). (14) 

The radius of convergence of this expansion is 2 because there is one singularity in the 5- 
plane at 5 = —2. [One might think (incorrectly) that as 5 ^ — 1 the integral for Z in Eq. (^ 
would become a representation of a Dirac delta function, which is singular. However, it is 
not generally true that when an integral representation for a function ceases to exist, the 



function exhibits a singularity [0.] Gat is analytic for Re(5 > —2 because the integration 



path in Eq. (|^) is an implicit function of 6. Indeed, as delta ranges through real values, the 
paths of integration of the two integrals in the first line of Eq. (|^ rotate in opposite directions 
||12|| . The path of integration of the first integral connects and oo in the complex-x plane 



along the straight line 



path 1 : argx = 37r/2 - 7r/(2 + (5). 



(15) 



The second integration path runs from to oo along 



path 2 : argx = 7r/(2 + 5) - 7r/2. 



(16) 



When 5 = (free field theory), path 1 connects — oo to and path 2 connects to oo 
along the real-x axis. For this case is real and there is no symmetry breaking. As 5 
increases, path 1 rotates anticlockwise and path 2 rotates clockwise. The two paths slope 
downward at 45° angles when 6 = 2 (0^ field theory). As 5 — > oo, both paths connect the 
origin to —too. Since the paths overlap and the integrations are in opposite senses, the two 
integrals in Eq. (0) cancel. [This cancellation becomes evident if we set 5 = oo in Eq. (|^).] 
However, for this theory the expectation values approach finite limits: 



The case 5 < is more interesting. As 6 decreases below 0, path 1 rotates clockwise 
and path 2 rotates anticlockwise. (The paths rotate infinitely fast as 5 ^ —2.) Whenever 



(5— >oo 



lim G 




(17) 



the two paths are horizontal, the resulting 02/(2fc+i) ^j-^ggry is real and there is no symmetry 



breaking. This happens for the special set of values 



6 



4k/{2k+l) {k = 0, 1, 2, 3,...). 



(18) 



For this case we have 




(19) 
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If the two paths overlap (this occurs when the paths are vertical and go up or down the 
imaginary axis), there is cancellation. This happens for a (j?-^^ theory when 

5 = {l-2k)/k (A; = 0, 1, 2, 3,...). (20) 

For this case we have 

lim G^ = (-.)^(-l)'=^[A;(iV+l)]!A!, (21) 

0— >(1— 2fe)/fe 

which reduces to Eq. ([T7|) for the special case = 0. 

Next, we examine theories in one- (^^mens^ona/ space-time (quantum mechanics). In gen- 
eral, the Hamiltonian 

H = pV2 + {ixf+^ (22) 

represents a broken-symmetry theory because H is not invariant under the parity operation 
P, whose effect is to make the replacements p —p and x —>■ —x. (While H is not symmetric 
under time reversal T, which makes the replacements p —p, x — x, and i —i, H is 
symmetric under VT.) We emphasize that there are both broken- and unbroken-symmetric 
phases in zero- and one-dimensional Euclidean space-time, but there is no transition between 
these different phases. Transitions can only occur when the dimension is two or more. 

As in the case of zero dimensions, there is a set of values of 6 [see Eq. (|18|)] corresponding 
to quantum-mechanical theories for which there is no broken symmetry. These special 
theories are described by the Hamiltonians 

= pV2 - a;2/(2fc+i) (A; = 0, 1, 2, 3,...), (23) 

where we have chosen = _x_ xhe Hamiltonians in Eq. (|23D have a deep connection 

with quasi-exactly solvable theories in quantum mechanics |]13[. To express our theories 
in a more familiar form we make the change of independent variable x = ^2^^+^ in the 
Schrodinger equation Hip{x) = Eip{x), followed by the multiplicative change of dependent 
variable ip = t^y to eliminate one-derivative terms. We obtain the Schrodinger equation 

7 



y"(t) = [(P + k)t-^ - 2{2k + l)H^\e + E)] y{t). (24) 

Equation represents a particle of zero energy in a rotationally symmetric quasi-exactly 
solvable potential having a centrifugal barrier. 

Finally, we use our new delta-expansion to calculate the one-point Green's function Gi 
in a symmetry-broken scalar quantum field theory in Z)-dimensional Euclidean space. The 
vacuum-persistence amplitude for the Lagrangian density in Eq. @) is 

Z = j p0e/^''"[-5(^'^)'-5-'</''+^'W)'+1. (25) 

Since the field has dimensions M^^~'^^^'^ and the coupling constant g has dimensions 
^2+6~SD/2^ where M is a mass, we re-express Eq. (^) using dimensionally explicit quantities 
by substituting g = ^M'^^^^^^^'^ . Expanding the exponent to first order in 5 gives 



X log (z(/)(t)M^-T) 1 +0(5^), (26) 



where 5o = | / d^x[{V(j)f + /i^^] and /i^ = + M^. 

To evaluate the one-point Green's function Gi we insert a factor of 0(0) into the integrand 
of Eq. (^) and divide by Z. By symmetry, the resulting expression vanishes at 5 = 0. To 
first order in 6 we have 

/ V(j) e-^«0(O) / d^t (f)\t) log (z0(t)Mi"f ) 

We now show that Gi is pure imaginary: We average the above expression with that 
obtained by setting to —0 at every lattice point. Provided that the integrals exist, we 
obtain the explicitly imaginary result 



r ■ ,,. JP0e-^»0(O)/ci^t0(t)|0(t)| 



It is not easy to evaluate this functional integral because it contains |0(t)|, so we employ 
a limiting process: lima^i/2 0^" = |0|- Before taking the limit we temporarily treat a as an 
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integer. Thus, this hmiting process involves continuing off the integers. [We cannot prove 
the vahdity of this technique, but we justify this procedure in part by comparing the final 
result with the first term in Eq. (0) for the case D = 0.] To evaluate this limit we introduce 
an external source J{x), and express Gi in terms of functional derivatives with respect to 
Jix): 

5 



where 



IT 



1+2q 



X lim 



Z[J] 



\SJ{t) 



(28) 



,/=o 



Z[J] = J p0e/'^''^[-^(^'^)'-5^'*'+^*] 



(29) 



and A(x, y), the free propagator in D-dimensional Euclidean space, is the Fourier transform 
of l/(p2 + /i2): 



A{x,y) 



-ip-{x-y)_ 



1 / 2tt\x — y\ 
2^ V Jl 



D 



K_^^D{fx\x -y\). 



(30) 



For integer a we use Eqs. (^) and (|30D to evaluate Eq. (|g) and obtain the final result 
for Gi: 

l+2a 



Gi = -i-SM^ lim / d^t 

4 a^l/2j 



S 



S 



6J{0) \6J{t) 



X 



d'^xd^yJ{x)J{y)A{x,y) 



a+l 



/(« + !)! 



.TT 2 (2a + 2)!A"(0,0)/rf^tA(0,t) 

—i—oM hm — 

4 a^i/2 2°+i(a + l)! 

-i(5My7rA(0,0)/2 j d^tA{0,t) 



?(5MV"'(/^V4vr)''/V^r(l - D/2)/2 



(31) 



where we have used / d t A(0,t) = /i and 



A(0, 0) = 12^-^2- ^7T-^/^T{l - D/2) {D < 2). (32) 
At D = 0, m = 0, and M = /i = a/2, Eq. (|31|) reduces to the first term in Eq. Note that 



as D increases past 2, Gi in Eq. ( pT]) becomes rea/ because the gamma function changes sign. 
For D < 2, A(0, 0) is finite and positive, but for > 2 it is a divergent integral regulated by 
continuing through complex dimension. Thus, for D > 2, the sign is a priori unpredictable. 
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